function z = find_omega_N1(g1)
global e_real tauL tauK betta alpha depr sigmaC sigmaL
kss=(((1/betta-1)/(1-tauK)+depr)/alpha/e_real^(1-alpha))^(1/(alpha-1));
css=kss^alpha*e_real^(1-alpha)-g1-depr*kss;
omega=(1-alpha)*(kss/e_real)^alpha*(1-tauL)*css^(-sigmaC)/(e_real)^(sigmaL);
z=(css^(-sigmaC)*css-omega*e_real^(sigmaL)*e_real)/(1-betta)-css^(-sigmaC)*kss*(1+(alpha*(e_real/kss)^(1-alpha)-depr)*(1-tauK));
